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ABSTRACT 

A ground-based polarimeter, PIQUE, operating at 90 GHz has set a new limit on the magnitude of 
any polarized anisotropy in the cosmic microwave background. The combination of the scan strategy 
and full width half maximum beam of 0?235 gives broad window functions with (£e) — 211^1% an( i 
(is) — 212^35 for the E- and B-mode window functions, respectively. A joint likelihood analysis yields 
simultaneous 95% confidence level flat band power limits of 14 and 13 /iK on the amplitudes of the E- 
and B-mode angular power spectra, respectively. Assuming no B-modes, a 95% confidence limit of 10 
fiK is placed on the amplitude of the E-mode angular power spectrum alone. 

Subject headings: cosmology: cosmic background radiation — cosmology: polarization — cosmology: 
observations 



1. INTRODUCTION 

Observations of the temperature anisotropies of the 
cosmic microwave background (CMB) already constrain 
cosmological models (e.g. Tegmark & Zaldarriaga 2000). 
In addition, a linearly polarized component of the CMB 
arises from any quadrupolar variation in the photons 
scattered from electrons at the last scattering surface 
(Rees 1968; Hu & White 1997). Polarization data will 
complement other CMB data (e.g., Kosowsky 1999). 
Current cosmological models predict that polarization 
anisotropies are 10-20 times smaller than temperature 
anisotropies. To date only upper limits on the CMB 
polarization exist, summarized in Staggs et al. (1999) 
and Subrahmanyan et al. (2000). The polarization 
fluctuations may be decomposed into gradient and curl 
parts (Kamionkowski et al. 1997), here designated E- 
and B-modes (Zaldarriaga & Seljak 1997). The complete 
formalism for analyzing polarization data developed 
by these authors, and extended in Zaldarriaga, 1998 
(henceforth, Z98), is applied to data here for the first time. 

2. INSTRUMENT 

The Princeton IQU Experiment 4 (PIQUE) comprises a 
single 90 GHz correlation polarimeter underilluminating a 
1.2 m off-axis parabola (Wollack et al. 1997). The beam 
on the sky is 0?235, and the instrument observes a single 
Stokes parameter in a ring of radius 1° around the NCP. 
The receiver is a heterodyne analog correlation polarimeter 
using W-band (84-100 GHz) HEMT amplifiers. A 
mechanical refrigerator cools the corrugated feed horn, the 
orthomode transducer (OMT), and the HEMT amplifiers 
to < 40 K. The telescope is fixed in elevation, but rotates 
in azimuth. Two large nested ground shields surround the 
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instrument; the inner shield rotates with the telescope. 

The correlation polarimeter (Krauss 1986) directly 
measures the polarized electric field, rather than detecting 
and then differencing two large intensity signals. The two 
input signals come from an OMT oriented so one arm is 
parallel to the azimuthal scan direction. The polarimeter 
output is proportional to the linear polarization for axes 
rotated 45° with respect to the inputs. The local 
oscillator signal in one polarimeter arm is phase-switched 
at 4 kHz, well above the 1 kHz 1// knee of the 
amplifiers (Wollack & Pospieszalski 1998). The 2-18 GHz 
intermediate frequency (IF) band is split into three 
frequency channels (SO, SI, S2). The total power in each 
of the IF arms is also detected, though with significantly 
less sensitivity since the phase-switching does not apply. 
The polarimeter characteristics are given in Table 1. 

3. OBSERVATIONS AND CALIBRATION. 

Between 2000 January 19 and 2000 April 2, PIQUE 
collected slightly over 800 hours of data from the roof of 
the physics building at Princeton, latitude 40?345, east 
longitude — 74?647. Every five seconds, the telescope 
alternated between two azimuth positions, ±0?93, at 
elevation 41?04, measuring =fQ (as defined by the IAU) 
for two regions separated by six hours in right ascension 
(RA) on the 6 = 89° ring. 

Since the polarimetry channels have a small (< 23 dB) 
sensitivity to total power, PIQUE differences data taken 
in the east position from data taken in the west postion 
five seconds later. Residual offsets in these "chopped" 
data, listed in Table 1, are small and stable (1// knee 
undetectable at 10 /iHz in clear weather). The offsets 
have been traced to polarized emissive pickup from the 
fixed ground screen, and are removed as described below. 
Since Q is measured in the west and — Q is measured in the 
east, the chopped data comprise sums of Stokes parameter 
Q. 

Constant elevation scans of Jupiter are used to 
determine pointing accuracy, map the beams, and 
calibrate the total power channels. The absolute errors in 
elevation and azimuth are 0?01 and 0?02, respectively. The 
beam FWHM are 0?233(9) (cross-elevation) and 0?240(14) 
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(elevation), in agreement with calculations and near- field 
data. For the purpose of the likelihood analysis, the 
beam dispersions are taken as a = 6*™ HM /V8 In 2 = 
0?1. Uncertainty in the brightness temperature of Jupiter 
contributes to a 10% error in the calibration of the 
total power channels (which are only used to correct for 
the slowly varying atmospheric opacity during the CMB 
polarization observations.) 

The polarimetry channels are calibrated by observations 
of the polarized emission (and reflection) from a large 
ambient-temperature aluminum flat. Nutating the flat 
about a vertical axis results in a peak-to-peak polarized 
signal of about 30 mK. These measurements are consistent 
with polarized cold load tests. The final calibration error 
is 10%, dominated by uncertainty in the surface resistivity 
(4.0 juO-cm) of the nutating flat. 

4. DATA REDUCTION 

Of the 807 hours of data, 200 hours were spent 
slewing the telescope, 123 hours were corrupted by various 
electromechanical failures, and 52 hours were isolated 
segments less than eight hours in duration which are 
not used. High levels of atmospheric noise contaminated 
102 of the remaining 432 hours, as determined by the 
large positive tails in the distributions of the polarization 
channels' correlation coefficients. In order to avoid biasing 
the data set, the primary atmosphere cut employs the 
quadrature chop data. A quadrature chop has the raw 
data from each azimuthal position (east and west) split in 
two to obtain (E\ + W\ — E 2 — W^jl, which contains no 
polarization signal. This cut removes 78 hours by requiring 
the absolute value of all three correlation coefficients for 
ten-minute sections of the quadrature data to lie below a 
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a Effcctive center frequencies and bandwidths, including gain slope 
and phase match degradations. 

^Raw sensitivity measured on a clear day with T a tm ~ 45 K, 
uncorrected for atmospheric opacity. The differencing strategy 
degrades the sensitivity by a factor of two. 

c Average chopped offset where the error gives the standard 
deviation of the mean of the eleven offsets calculated during the 
season; see text. The typical statistical error in calculation of a 
single one of the 11 offsets is 20 /iK. 

d The average error per bin for the final data binned into 
twenty- four 0? 26- wide RA bins. Note that the analysis uses data 
in 144 bins. 



Note. — All temperatures are in thermodynamic units. 
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Fig. 1. — Binned data in thermodynamic units for each 
frequency channel, versus RA. For this plot, 24 approximately 
beam-sized bins are shown, though the analysis uses 144 bins. 
The theoretical point source response to a 250 mjy polarized 
source is indicated in the top panel. 



threshold, r/ q = 0.32. Varying r) q by ±5% does not affect 
the final results statistically. The secondary atmosphere 
cut culls extreme values (> 4tr) of the CMB correlation 
coefficients and removes an additional 24 hours of data. 

The final data cut makes use of the six hour null test. 
The six hour null test takes advantage of the fact that 
certain linear combinations of data separated by six hours 
in RA should have zero signal: d t — d t +Q+d t +i2 — d t +is = 0, 
where t is measured in hours. Those days that fail this null 
test (p{x 2 W) < 0.05) are removed. This final cut removes 
an additional 82 hours of data, leaving 248 hours of CMB 
data for the final analysis. Note that while this cut only 
slightly increases the final limits, it significantly improves 
the H1-H2 null test described below. 

The remaining time series is binned into w 2100 
10-minutc-wide bins. The means, standard deviations and 
intcr-channcl correlations are calculated for each bin. The 
attenuation due to atmospheric opacity for each bin is 
estimated from a total power channel and the appropriate 
(w 15%) correction is applied. The data are divided into 
eleven periods separated from one another by stretches of 
weather severe enough to require covering the instrument. 
The periods have varying lengths t,, where 16 hr < t, < 
110 hr. An offset for each channel is calculated and 
removed from each period. Other divisions of data, in 
which as few as two or as many as 24 offsets are removed, 
do not change the limit by more than 1 ^K. The weighted 
means, standard deviations, and covariances are calculated 
for 144 10-minute RA bins. These data are rebinned for 
presentation in Figure 1. 

Table 2 presents the results of a series of null tests. The 
null data sets include the quadrature data, data differenced 
between the first and second halves of the observing run 
(H1-H2), data differences between signal channels, and 
the six-hour null data. The \ 2 distribution of the null 
tests is consistent with noise. Similarly, no contamination 



is found when the data are binned in Sun-centered or 
Moon-centered coordinates. 

5. DATA ANALYSIS 

The Bayesian method with uniform prior is used to 
establish limits on the amplitude of the polarization 
anisotropy and to verify null tests of the data. The data 
vector d has 3N entries, where N = 144. The data 
covariance matrix Co consists of nine N x N diagonal 
submatrices, since correlations between 10-minute-wide 
bins are 
small enough to neglect (<C 1% of the variances). Six 
of the submatrices account for inter-channel correlations 
due to HEMT-induced correlated gain fluctuations and 
correlations from atmospheric fluctuations weakly coupled 
into the polarimeter channels. The average correlation 
coefficients between S0-S1, S0-S2, and S1-S2 are 0.02, 0.08 
and -0.01, respectively. The theory covariance matrix Ct 
may be expressed in terms of the E- and B-mode angular 
power spectra Cf and Cf by 
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(2£ + l)[Cf + CfWi^a)], 



where fiij are position vectors, fcj = cos 1 (hi ■ hj) is the 
lag, and Wi t e and W2,t are associated window functions. 
The Q indicates that PIQUE measures sums of Q's. Z98 
presents an elegant way to describe the window functions 
for a ring observation strategy, here modified to include 
the PIQUE differencing: 



where a = tt/2 is the chop amplitude, and = 7r/72 
is the bin size. The window functions are depicted in 
Figure 2. The likelihood is C oc exp(-d T M~ 1 d/2)/ y/\M~\, 
where, M = C T + C D . 

The likelihood analysis proceeds by first considering flat 
angular spectra, such that £(£ + \)Cf /2n — Tj, where 
X = E,B. Since the amplitude of Cf is predicted 
to be much smaller than that of Cf (see Figure 2), 
Table 3 presents the limit on Te under the assumption 
T B = 0. Figure 3 shows C(T E ,0) for the CMB data 

and two null data sets. The limit Te is found by 
integrating £(T E , 0); the result is w 30% higher if £(T§, 0) 

is integrated. Next, joint upper limits (T e ,Tb) are 
determined by finding the constant contour of C(T e ,Tb) 
enclosing 95% of the volume of C Finally, a gaussian 
autocorrelation function (GACF) analysis is performed. 
The motivation to revive the GACF is that it better 
describes the gross features of the predicted E-mode 
spectrum than the flat model; see Figure 2. For the GACF 
with characteristic scale £ c , the power spectra are given by 
£{£ + l)Cf/2n = A 2 w 2 cxp(-ii 2 /2), where X = E ,B Q , 
u w £/£ c (Bond 1995), and £ p = V2£ c is the peak of the 
GACF. Only the case Bo = is considered. The likelihood 
£(E 0} 0) is integrated to place 95% confidence level upper 
limits on E for several values of £ p as shown in Figure 2. 
Table 3 summarizes the limits obtained for the CMB data 
and two null tests. The limits in Table 3 do not include 
calibration errors (10%). Note that for PIQUE's broad 
window functions, the 7% beam errors only add 2% errors 
in quadrature with the calibration errors. 



W {lt2 }A4>ij) = 2 E F {i,2},em(0)BL cos(m^), (2) 

m— 1 

where the factor of 2 in front accounts for the negative m 
values, and the m = term is excluded from the window 
functions because offsets are removed from the data. The 
F{i,2},tm are given in Z98 in terms of associated Legendrc 
polynomials, — 1° is the radius of the ring, and 



GACF Limits on E„ _ 
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Results of \ 2 consistency tests. 
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The data sets are described in the text. 
^The number of degrees of freedom for each data set. 
"The probability of exceeding the \ 2 f° r a given frequency channel. 
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Fig. 2. — Zero-lag window functions for E- and B-modes 
(solid and dashed lines) are plotted in the lower panel. The 
middle panel shows the limits on Tg and Tb described in 
the text, plotted over E-mode predictions (dash-dot) from the 
Qtot = 1 model of Jaffe et al. (2000) and B-mode (solid, 
assuming T/S = 1). The smooth dashed curve is a GACF 
with E = 5.6 nK and 4 = 687. To illustrate PIQUE's ability 
to limit the fluctuations where they are thought to be sizeable, 
the top panel shows the 95% confidence limits on the amplitude 
of the GACF, Eo, where the abscissa is l v = \[2£ c . 
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Normalized likelihoods for the cases when Tb = 0. 

6. DISCUSSION 

The formalism of Z98 enables us to place limits on both 
the E- and B-mode angular power spectra, allowing direct 
comparison of these results with theoretical expectations. 
Due to PIQUE's broad reach in f-space, these results limit 
the size of the predicted acoustic peaks in the polarization, 
rather than the damped tails at high and low I. PIQUE's 
limits on Te and Tb are higher than expected from, e.g., 
the Q tot = 1 model of Jaffe ct al. (2000). The limits 
are also higher than expected foreground levels. Polarized 
dust emission anisotropy, extrapolated from the IRAS 
100 /xm map (Whcclock et al. 1994), should contribute 
< 0.5 ^K, while polarized synchrotron, extrapolated from 
Brouw & Spoelstra (1976) or from Haslam et al (1982), 
should be less than 0.4 ^.K. Spinning dust emission is even 
smaller (Draine & Lazarian 1999). 

Since this paper represents the first effort to put limits 
on Te and Tb directly, it is not possible to quantitatively 
compare these results to published results from previous 
experiments. Direct comparison would require other 
authors to reanalyze their data for sensitivity to E- and 
B-modes rather than to Q and U. Qualitatively, we note 
that the limit Te < 10 /jK is smaller by more than a 
factor of two than the best limits previously published 
for I < 3000 (Torbet et al. 1999; Netterfield 1995). 
Further, PIQUE has sensitivity to higher I than these two 
previous experiments, so this work represents a significant 
improvement in the constraints on the polarization of the 
CMB. 

We are indebted to Matias Zaldarriaga for insights on 

Table 3 
95% Confidence Limits 
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a The data sets are described in the text. 
b The limit Te is found assuming Tg = 0. 

c The limits Te and Tb are determined simultaneously by finding 
the contour of constant likelihood enclosing 95% of the volume. 



The limit on the amplitude of the GACF, Eq, at i v 
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